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Abstract 

The vacuum energies corresponding to massive Dirac fields with the 
boundary conditions of the MIT bag model are obtained. The calcula- 
tions are done with the fields occupying the regions inside and outside the 
bag, separately. The renormalization procedure for each of the situations 
is studied in detail, in particular the differences occurring with respect to 
the case when the field extends over the whole space. The final result 
contains several constants undergoing renormalization, which can be de- 
termined only experimentally. The non-trivial finite parts which appear 
in the massive case are found exactly, providing a precise determination 
of the complete, renormalized zero-point energy for the first time, in the 
fermionic case. The vacuum energy behaves like inverse powers of the mass 
for large masses. 
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1 I. Introduction 



The first modern calculation of the vacuum energy density of a quantum field 
in the presence of boundaries is almost fifty years old. As is well known, it is 
due to H. Casimir 0. Its first measurable consequence was the attraction in the 
electromagnetic vacuum of two neutral, infinitely conducting plates (thereafter 
called Casimir effect, see for instance Ref. 0]). There had been, before that 
time, other calculations and explanations of the attraction of two neutral bod- 
ies, understanding them as mere van der Waals effects, ||, and also a paper by 
Casimir and Polder [Q where the finiteness of the velocity of light had been taken 
into account. Casimir's paper |I] was the first one where an absolutely modern 
quantum field theoretical calculation was performed, using the concept of zero- 
point energy (whose physical relevance was somehow unclear at that time). The 
treatment of the divergences resulting from the infinitely many degrees of free- 
dom has been the most difficult point. Since then, calculations of the vacuum 
energy have attracted the interest of many scientists because it turns out that, in 
different contexts, the inclusion of quantum fluctuations about semi-classical con- 
figurations is essential. On the other hand, spherically symmetrical situations are 
very important for practical applications. The calculations involved are certainly 
much more complicated than in the case of systems with plane boundaries. 

Historically the first far reaching ideas involving vacuum energies in the case 
of spherical configurations also originated with Casimir. He proposed that the 
force stabilizing a classical electron model arises from the zero-point energy of the 
electromagnetic field within and without a perfectly conducting spherical shell 
||. Having found an attractive force between parallel plates due to the vacuum 
energy [l] , the hope was that the same would occur for the spherically symmetric 
situation. Unfortunately, as Boyer |J first showed, for this geometry the stress is 
repulsive |7], ||. Nowadays it is known that the Casimir energy depends strongly 
on the geometry of the space-time and on the boundary conditions imposed. This 
is a very active field of research (see, for instance, H UPl ). 



More recently the zero-point energy has received considerable attention in 



the context of the bag model |TTJ-[|T5| and chiral bag model |16|-||22||. In these 
systems, quarks and gluons are free inside the bag, but are absolutely confined to 
it, being unable to cross the boundary surface. This is imposed, mathematically, 
by appropriate boundary conditions. The sum of the mesonic, valence quark 
and vacuum quark contributions to the baryonic number have been found to be 
independent of the bag radius and pion field strength, being the vacuum quark 
contributions — which are analogues of the Casimir effect in QED — essential in 
the calculation of baryonic observables. The issue of regularization in this model 
is certainly non-trivial. Under specific circumstances, different regularization 
procedures can yield different results and real physical problems arise in the 
calculation of quark vaccum contributions to some barionic observables, as the 
energy itself. 
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The Casimir energy for a spherical capacitor is usually given by the sum of two 
terms: an internal and an external one. Volume divergences cancel in each of the 
two regions separately, but surface and curvature divergences survive in each part 
and only cancel when interior and exterior contributions are added up, leaving 
then constant terms. In chiral bag models an additional divergence -logarithmic 
in the plasma frequency — appears. Having done these considerations, one should 
observe, however, that in the bag model quarks are supposed to exist in the 
interior of the bag only and, therefore, there is no clear way of eliminating these 
surface, curvature and logarithmic divergences. Some authors (see have 
suggested a revision of the Cheshire cat principles in the sense of including free 
quarks as correct high-energy degrees of freedom for the bag exterior, which 
would appear with a smooth transition at energies of the order of some GeV. 
In this way divergences would cancel quite naturally. However, in our opinion 
this deconfining transition is even a harder issue than those of regularization and 
renormalization. It will not be dealt with here. 

For massless fermions the zero-point energy was considered some time ago in 



Refs. p3| , p4| , finite temperature effects were taken into account in [p~|, |2~I| . The 



massless fermionic field inside and outside the spherical bag was analysed in [g4fl . 
In the last mutual cancellation of the divergences of the inner and outer 

spaces occurs. As a result, finite zero-point energies are found. However, when 
considering only the inner space, divergences arise and it is necessary to introduce 
contact terms and perform a renormalization of their coupling. Results for the 
massive fermionic fields contain new ultraviolet-divergent terms in addition to 



those occurring in the massless has been discussed in ||25|| . Further 

considerations, especially on the renormalization procedure — necessary in order 
to carry out these calculations — and also on its precise physical interpretation, 



can be found in [26 



In most of the papers mentioned above a Green's function approach has been 
used in order to calculate the zero-point energy. An exception is Ref. fl2"6"| , where, 
in the general setting of an ultrastatic spacetime with or without boundaries, a 
systematic procedure which makes use of zeta function regularization was de- 
veloped. In this approach, a knowledge of the zeta function of the operator 
associated with the field equation together with (eventually) some appropriate 
boundary conditions is needed. Recently, a detailed description of how to ob- 
tain the zeta function for a massive scalar field inside a ball satisfying Dirichlet 
or Robin boundary conditions has been given elsewhere by the authors of the 
present work |27|, |28|. An analytical continuation to the whole complex plane 
was obtained there and subsequently applied to the computation of an arbitrary 



number of heat-kernel coefficients. In ensuing papers [29, p(J the functional de- 



terminant was considered too and, furthermore, the method has been also applied 
to spinors |3T| and p-forms ^ |55| . All the above considerations are purely 
analytical and quite precise. In order to obtain explicit values for the Casimir 
energy, however, a numerical evaluation of an integral was necessary. This has 
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been done in different cases, in particular for the massless scalar field and the 
electromagnetic field |37|], partly reobtaining previous results. 

To finish this description of recent previous work, let us mention that in 



Ref. |38| we have investigated the case of a massive scalar field in the bag. 
We have discussed there how, for the case of a massive field — already for a 
scalar one — non-trivial finite parts which depend on an adimensional variable 
involving the mass are present, that need to be properly renormalized, in order 
to get the corresponding zero-point energy. In the present paper we shall extend 
our analysis to the case of Dirac fields, generalizing in this way our considerations 
to a situation that approaches very much the conditions of a realistic MIT bag 
model. 

The organization of the paper is as follows. We shall rely on our previous work 
(dealing with the bosonic case) for a precise description of our method — which 
was given there in full detail [|3S| — as well as for the particular formulas that are 
needed in the subsequent study of the zeta function of the problem we consider 
here. We feel that to repeat all this here would not be justified. Consequently, 
in Sect. 2 we will proceed already with the specific description of the model for 
the case of Dirac fields inside the bag with boundary conditions corresponding 
to the MIT bag. Starting from the Dirac equation and imposing the boundary 
conditions we will derive an eigenvalue equation in terms of Bessel functions. 
This will be the basic equation to solve, what we shall do in the same section for 
the interior of the bag. In section 3 we will explain the renormalization scheme 
used for the model. Sect. 4 contains the analogous treatment for the region 
exterior to the bag and for the whole space. Adding up the interior and the 
exterior contributions, we will see how the divergences cancel among themselves, 
as well as the influence of this cancellation on the compulsory renormalization 
process. It turns out that important differences with respect to the non-fermionic 
case appear concerning this issue, although we shall argue that, in the end, they 
will not affect substantially the interpretation of the physical results. Sect. 5 is 
devoted to conclusions. The appendixs contain some hints and technical details 
that have been used in the derivation of the zeta function (App. A) and a full 
list of the constituents that build up the subtraction terms in the decomposition 
of the zeta function, an essential (although rather technical) step in our method 
(App. B). 



2 Fermions inside the bag 

The first task is to derive the energy eigenvalue equations for a Dirac spinor 
subject to the MIT bag boundary conditions. The setting we consider first is the 
Dirac spinor inside a spherically symmetric bag confined to it by the appropriate 
boundary conditions. The coordinates we use are just the spherical ones, r, 6, tp, 
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which best adapt to the form of the bag. Thus, we must solve the equation: 



H(P n (r)=E n n (r), (2.1) 

H being the Hamiltonian, 

H = -«7 V<9 Q + 7°m, (2.2) 
with the boundary conditions 



1 + i -7 



0. (2.3) 



These boundary conditions guarantee that no quark current is lost through the 
boundary. 

The separation to be carried out in the eigenvalue equation Q2.ip is rather stan- 
dard and will not be given here in detail. Let J be the total angular momentum 
operator and K the spin projection operator. Then there exists a simultane- 
ous set of eigenvectors of H, J 2 , J 3 , K and the parity P. The eigenfunctions for 
positive eigenvalues k = j + 1/2 of K read 



\ - yi^(u,r)CV m (J) j ' 



(2.4) 



whereas, for k = —(j + 1/2), one finds 

(2.5) 



Here = \/ E 2 — m 2 , A is a normalization constant and Qji m (f/r) are the well 
known spinor harmonics. In order to obtain eigenfunctions of the parity operator 
we must set I' = I — 1 in ( |2.4| ) and I' — I + 1 in (|2.5|). In both cases, j = 
1/2, 3/2, oo, and the eigenvalues are degenerate in m = —j, +j. 

Imposing the boundary conditions ( |2.3|) on the solutions and Q2.5|) , re- 
spectively, one easily finds the corresponding implicit eigenvalue equation. For 
k > 0, it reads 



V E — m 
and for k < 0, on its turn, 



/: m J j+ i{ujR) + Jj{uoR) = 0, (2.6) 



JjM) - \/§T^ Ji+l(a;il!) = °' (2,7) 
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Regretfully, it is not possible to find an explicit solution of equations ( |2.6|) and 
([2.7]) . But as we have shown in our previous paper for the case of the scalar 
field — and will describe below for the spinor field — the information displayed in 
( |2.6|) and ( f2T7"D is already enough for the calculation of the ground state energy 
for massive spinors in the bag. 

The regularization of this ground state energy will be performed by using the 
zeta function method. In short, we consider 

Eo(s) = -^EK^/A Res>s = 2 



-A 



(int) i 



1/2)// 



2.s 



(2i 



and later analytically continue to the value s = in the complex plane. Here so 
is the abscissa of convergence of the series, /i the usual mass parameter and 



k 



(2.9) 



The power of the method lies in the well defined prescriptions and procedures 
that we have at our hand to analytically continue the series to the rest of the 
complex s-plane, even when the spectrum is not known explicitly (as will in 
fact be the case). These procedures have been developed and described in great 
I, so that we can be brief here. 



detail in [27 



The zeta function in the interior space is given by 



c M) (s) 



oo „ 

2 E (2J + 1)/ 

J=l/2,3/2,... J 1 



dk 
2tt7 



9 i 



(k 2 + m 2 Y 



2m 



Jj(kR) - J- +1 (kR) + —J j (kR)J j+1 (kR) 



(2.10) 



Here the factor of 2 results from taking into account particles and antiparticles. 
Using the method — ordinarily employed in this situation — of deforming the 
contour which originally encloses the singular points on the real axis, until it 
covers the imaginary axis, after simple manipulations we obtain the following 
equivalent expression for Q mt >: 



c 



(int) i 



2 sin 7ts 



E (2J + 1) / dz 

■ ,7T!,, n JmR i 



TV 

oz 

( m -^) / J(^') / /(^*) 




— m 



(2.11) 



2mR\ ,, 2/ , 



z 2 j 



5 



As is usual, we will now split the zeta function into two parts: 



^ int \s) = Z N (s) + E Ms) 



N 



(2.12) 



namely a regular one, Zn, and a remainder that contains the contributions of 
the N first terms of the Bessel functions I v {k) as u, k — > oo with vjk fixed [39]. 
The number N of terms that have to be subtracted is in general the minimal 
one necessary in order to absorb all possible divergent contributions into the 
groundstate energy, Eq. (|2.8| ). In our case, N = 3. This is a general procedure, 
commonly applied in order to deal with such kind of divergences. We get 



Z s (s) 



sinus J^, . f°° , \/zjy 



m 



J z a z Zj J zj R J 



In 



e^l + z^l-t) 



* D k {mR,t) 

2-^ Ak 

k=l J 



(2.13) 



where 77 = \/l + z 2 + ln[2;/(l + yl + z 2 )] and t = 1/vl + z 2 . After renaming 
mi? = x, the relevant polynomials are given by 



Di(t) = - + (x-l/4)f 



t 6 t 5 



.r 



— - — + + 1/8 r + -- + 1/8 r - 



.r 



2^,2 



179 1 9 3t s 



576 



23 7s 



64 8 



9 



+ — + — — + — )f + (x - 1/2) t b +[ — -- + — ]t 



320 4 



+ 



x 



X 



X 



+ 1/8 - £ W 4 + ~ + 



8 192 



+ 



x 



t 6 



(2.14) 



The asymptotic contributions A^s), i = —1, ...,3, are defined as 

2 ^- s VT+x^-i 



A (s) 
Ms) 



6 sm{7rs) 

71 .=1/2 

4 sintVs^ x 



7T 



E j(j + i/2) / 

=1/2 J " 

00 . 

E (7 + 1/2) / 



rnR/j 



■II 

R 



m 



x 



3=1/2 



inR/j 



— I — m I — In 



d v / TTx2(l -t) 



^E0 + i/2) 
77 ,=1/2 



mJi/j \\ R J 



m 



dx 

d Djjt) 
dx j i 



x- 



(2.15) 
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The small mass expansion can be conveniently represented as 



A-i(s) = 



d2s oo 

it \ - 



k\ 



{mR) 



k + s 



x 



2C(2fc + 2s-2,^) + C(2A; + 2s-l,^ 



A (s) 



R 



2s 



E 



-v. 



(mR) 



2k 



T(k + s + \) 



A(s) = 



2v / ^r(s)^o k\ k + s 

x (2C(2fc + 2s-l,^) + C(2A; + 2s,^)) 

2R 2s oo 

"r5J 



fc=0 



fc! 



x 



* T(k + s + ^) 

/ , X i,a 
a=0 



r(s±i) 



2C(2A; + 2s + i - 1, -) + ((2k + 2s + ^) 



(2.16) 



In this expression, the x ita are the coefficients of the expansion of the functions 
Di(t), i.e., 



2i 



Di(t) = 5> ii0 £ 

a=0 



a+i 



(2.17) 



Note that here we encounter the same problem that occurred already in the scalar 
case. One needs a representation that is useful and valid for an (in principle) 
arbitrary value of m. To this end one can actually proceed in different ways, 
casting the final result in terms of convergent series or integrals. Our leitmotiv 
will be the following: we will always try to express the final result in terms of 
the formula which is more appropriate for practical evaluation (e.g., numerical, 
in general). This means that, sometimes, instead of having the closed convergent 
sums that could be universally used in the scalar case, rapidly converging integrals 
better suited for numerical analysis will be here preferred. 

With this aim, we note that after performing the z-integration the A^s), for 
% > 1, can be written in the following form, 



Ms) = 



\nr 2s ^ 



x i>a T(s + + a)/2) 



x 



T(s) £Z(mR)»* m + a)/2) 
f(s; 1 + a; (i + a)/2) + a; (i + a)/2) 



(2.18) 



with the definition 

f(s;a; b) 



E » a i + 

1^=1/2,3/2,... \ 



— V 

mR) 



-s-b 



(2.19) 
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The remaining thing to do in the present case is to calculate the f(s; a; b) for 
the relevant values at s — —1/2. This is a systematic calculation that will be 
sketched in App. A. Let us here mention only that an essential step is to use the 
simple recurrence: 

f(s; a; b) = (mR) 2 [f(s; a - 2; b - 1) - f(s; a - 2; b)} . (2.20) 

In App. B we give the whole list of starting terms that, in addition to the 
recurrence formula ( |2.20| ), are strictly necessary for obtaining explicitly all the 



AAs) needed in our calculation. 



3 Discussion of the renormalization 

For the discussion of the renormalization let us look for the divergent terms in the 
groundstate energy. By construction they are all contained in the contributions 
Ai(s). Having their explicit form at hand (see Apps. A and B) they can be given 
quickly. In particular, we have, for the interior part: 

res A ( lr (-1/2) = + + 



res 4" 4) (-l/2) = 

V ' 1 2vr 24ttR 

res A\ mt) (-l/2) = + + 



6tt 12tt 4807ri? 
m 2 R 1 



res 4 in<) (-l/2) = - m ( £ + ) + — — + 



7i 12tt 12vr 48?^ 
m 2 R /l 1 \ 1 



4 V8 2tx) 128R 24nR 

Aim), /rA 2m 3 R 2 j D /l 2 

res A\ —1/2 = h m it - H 

3 V 1 ' 3vr V4 3tt 

1 7 \ 1 97 
+m I - + 



8 20tt/ 128i? 100807ri? 



and, clS db result, 



bSirR 157T d7T 07T 07T 

These terms form the minimal set of counterterms necessary in order to renor- 
malize our theory. 

In the scalar case one had the peculiar situation that there were no divergent 
contributions of the form ~ m 3 , m in the zeta function description |38[] — although 
in other regularizations they indeed appear . So in principle one had the choice 
of renormalizing the associated couplings. Contrarily, for spinors, as seen in ( |3.1| ), 
the coupling constants of all terms appearing have to be renormalized. We are 
led into a physical system consisting of two parts: 
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1. A classical system consisting of a spherical surface ('bag') with radius R. 
Its energy reads: 

E class = pV + crS + FR + k + -, (3.2) 

ri 

where V = f tt-R 3 and S = 4nR 2 are the volume and the surface of the bag, 
respectively. This energy is determined by the parameters: p pressure, a 
surface tension, and F, k, and h which do not have special names. 

2. A spinor quantum field ip(x) obeying the Dirac equation and the MIT 
boundary conditions Q2.3|) on the surface. The quantum field has a ground 
state energy given by E , Eq. (|2.8|). 



Thus, the complete energy of the physical system is 

E = E dass + E Q (3.3) 

and in this context the renormalization can be achieved by shifting the parameters 
in E c i ass by an amount which cancels the divergent contributions. 

First we perform a kind of minimal subtraction, where only the divergent 
contribution is eliminated, 

m 4 1 m 3 1 



p — >• p - ttt^t- cr -> a 



167T 2 S 247T 2 S 

F - F+^- k^k-^ 1 - (3.4) 
6tt s 30n s 

1267TS 

The quantities a = {p, a, F, k, h} are a set of free parameters of the theory to be 
determined experimentally. In principle we are free to perform finite renormal- 
izations at our choice of all the parameters. 

We find natural to perform two further renormalizations. First it is possible to 
determine the asymptotic behavior of the A% for m — > oo using the results of Apps. 
A and B. The finite pieces not vanishing in the limit m —>■ oo are all of the same 
type appearing in the classical energy. Our first finite renormalization is such 
that those pieces are cancelled. As a result, only the "quantum contributions" 
are finally included, because, physically, a quantum field of infinite mass is not 
expected to fluctuate. The resulting A» will be called Af en \ 

Concerning Z3 we have not been able to determine analytically its complete 
non- vanishing behavior for m — >• 00. Instead, for the numerical analysis, as shown 
in Fig. 1, we have constructed a numerical fit of Z% by a polynomial of the form 

4 

P{m) = ^Cjjn 1 , 

i=0 
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and then subtracted this polynomial from Z3. As explained above, this is nothing 
else than an ulterior finite renormalization. The result will be denoted by Z^ ren \ 
Summing up, we can write the complete energy as 

E = E dass + Et n \ (3.5) 

where E c i ass is defined as in (|3.2|) with the renormalized parameters a and Eq^ = 

Figure 1 shows the numerical analysis of the energy Eq^ of the system for 
this specific choice of renormalization. The energy exhibits a clear minimum 
corresponding to a stability bag radius. 



4 Exterior of the bag and a model for the whole 
space 

The analysis of the region exterior to the bag is quite similar to the one carried 
out for the interior region. Only some specific differences appear both in the 
formulas and in the results. The expression for the zeta function in the exterior 
region is essentially the same as the one corresponding to the interior, but for the 
replacement of the Bessel Ij functions with Bessel Kj functions, namely 

C^\s) = £ (2j + 1) r . dz {{zj/Rf - m 2 



.7=1/2,3/2,..- JmR H 



9 1 

x 7T ln 

oz 



z 2j U}{*5) - Aj. ,(:./) + ^Kj(zj)K j+1 (zj) 



The splitting of the zeta function has also the same aspect as for the interior 
region. We have, in particular 

A { *?\s) = A { i?\s), 



4 sinfTTs) 



7T 



oo poo 

ECj + 1/2)/ , dz {{zj/Rf -m 2 ) 



3=1/2 
X ^~ ln 

oz 



l + t 



and the polynomials that replace the Di(t) above are here (x = mR) 
Dt(t) 



t t 3 

— h xt , 

4 12' 



— -x 2 t 2 t 3 xt 3 t A xt 4 t 5 t 6 

D 2 (t) = 1 1 1 , 

w 2 828288 
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— , , -5t 3 xt 3 

D 3 (t) = + 

v ; 192 8 

xt 5 x 2 t 5 t e R 23 1 7 7xt Y 3t s 179^ 



X 3 t 3 


t 4 




x t 4 x 2 1 4 


9t 5 




+ 8 


+ 




320 


3 


2 ' 2 


t 6 






23t 7 7xf 


3t 8 




-xt 6 


+ 




+ X 


" 2 




64 ' 8 



4 2 2 64 8 8 576 

As for the functions A\ (s), one obtains the same expressions as before, with 
the replacement of the polynomials D^t) with the corresponding polynomials 

A(t). 

In principle, the same procedure as before can be applied now in order to 
get an analytical expression for the whole energy of the exterior space. Instead, 
we want to restrict ourselves here to the specific changes that show up when 
discussing the renormalization. For doing that, we only have to consider the pole 
of the different A^\ xt ^ . In particular, we have for the residua 

-^(-1/2) = -^ + ^ + i ^ = res^)(- V2 ) 

resAtH-m = ^r + ^b? = - resA o m4) (-i/2) 

res 4"° (-1/2) = _!^!_!^ + m 1 



7i 12vr 12tt 487T.R 

,(ert) , , s m 2 R (\ 1 \ 1 I 

res 4 (-1/2) = ^ + m 



4 V8 2vr/ 128i? 247T.R 



Aext), , x 2m 3 R 2 9 (\ 2 \ 



1 7 \ 1 97 
-m | - 



20tt/ 128i? 10080tt.R 
This yields for the residue of the whole zeta function at the exterior region: 

Res c^(-V2) = - — + ^ (41) 

Thus the minimal set of counterterms necessary in order to renormalize the theory 
in the exterior of the bag is identical to the one in the interior of the bag. The 
classical system is again described by Eq. (|3.2|) . 

The opposite sign of the coefficients in the divergences ( |3.1| ) and ( |4.1| ) corre- 
sponding to the odd powers of R can be easily explained by means of differential 
geometrical arguments, just observing that the curvature of the surface of the 
bag has opposite sign when looked at from the exterior and from the interior of 
the bag. 

Contrary to the scalar case, the divergences from the two sides do not anihilate 
when adding up the two contributions. In fact, for the zeta function corresponding 
to the whole space (internal and external to the bag) we obtain: 

Res C(-l/2) = Res C^-l/2) + Res C (e ^(-l/2) = -— , (4.2) 

157T 37T 
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therefore, the two free parameters a and k remain even if the whole space is 
considered. The only exception is the massless field where for this reason the 
issue of renormalization is much simpler than in our case. 



5 Conclusions 

In this paper we have studied in considerable detail a quantum field theoreti- 
cal system of a Dirac field with boundary conditions corresponding to the MIT 
bag model. This is the most natural continuation — in the direction towards ap- 
proaching truly realistic physical systems — of previous work where only scalar 



fields were treated |38[. We have seen that the consideration of a fermionic field 
carries along a number of additional difficulties, mainly in relation with the phi- 
losophy of the renormalization process. Up to this point, the application of our 
techniques can be carried out essentially in the same way as for the scalar case. 
Starting from the Dirac equation and imposing the boundary conditions we have 
derived an eigenvalue equation in terms of Bessel functions. This basic expres- 
sion has then been solved, implicitly, in the regions interior and exterior to the 
bag surface, by using contour integration. This has yielded the corresponding 
zeta function in each of the two domains. Extraction of the singular part of the 
zeta function has been also done exactly. However, adding up the contributions 
of the two parts, not all divergences cancel among themselves (as was the case 
for a scalar field), what theoretically influences the playground of the ulterior 
renormalization process. 

In the end, a detailed consideration shows that, after renormalization, two 
dimensionless parameters remain whose values cannot be fixed theoretically, but 
have to be numerically adjusted by direct comparison with the physical system 
described by the model. In this, we must confess, we are still a bit far from our 
final goal. In the sense that, as it stands, our model cannnot be considered yet to 
describe a realistic physical situation. This must be left to future work, given the 
complexity of the proposal. In any case, we should like to point out the rigour 
and strict systematicity of the approach we have used here, and also its relative 
simplicty, if we compare it with other methods of similar strength and generality. 

As noticeable results of our analysis we would like to mention, that the Casimir 
energy may have a clear minimum associated with a stable bag radius (see Fig. 1). 
Comparing this behaviour with the one corresponding to the scalar field, where 
a maximum occured, this difference can clearly be traced back to the anticom- 
muting nature of the spinor fields, which shows up as a sign in the definition of 
the groundstate energy. 

Another interesting observation is that, contrary to the case of parallel plates, 
the behaviour of the Casimir energy for large values of mR is not exponentially 
damped. Instead, as is clearly observed from the representations of the Ai(s) 
given in the Apps. A and B, we find a behavior in inverse powers of the mass. 
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This is directly connected with the nonvanishing of the extrinsic curvature at the 
bag. 

Possible continuation of our approach is in the direction of finite temperature 

n. a 



and finite densities as considered already for massless fermions in [15, 21 



natural question to ask concerns the possible appearance of a first order phase 
transition from an hadronic bag to a deconfined quark-gluon plasma within our 
framework. This is left for future work. 
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A Appendix: Explicit representations for the 
asymptotic contributions inside the bag 



The essential formulas for the basic series f(s; a; b), eq. ( 2.19|) , in the calculation 
are the following: 



1^=1/2,3/2,. 



+ (-l)"2 / du — 1 

o 
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Using partial integrations one can get representations valid for values of s needed 
for the Ai(s). One gets, for example, the following expansions around s = —1/2: 
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2 1 + e 2 ™ 



: + 0(s + l/2) 



showing clearly that one can obtain quickly convergent integrals, respectively 
expressions for the effective numerical evaluation of the involved sums. 

All the particular values that are necessary to give the Ai(s), % = 1,2,3, 



explicitly (in addition to the recurrence ( |2.20|) ) are listed in App. B. 

The first two leading asymptotic contributions, A-i and A have to be treated 
in a slightly different way, as has been explained in detail in B3] . For completeness 
we give the final results 
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(A.3) 
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(A.4) 



This completes the description of our procedures to obtain numerically evaluable 
representations for all the A^s, needed for the calculation of the Casimir energy 
of the spinor inside the bag. 
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B Appendix: Full list of constituent terms /(a; b) 
to be used in addition to the recurrence for- 
mula 

To simplify the expressions, we shall here use x for mR. In addition to the above 
recurrence, in order to determine all the A^s) explicitly one needs the following 
f(a; 6)'s [we shall use the notation f(a; b) = /(— 1/2; a; b)]: 
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With them, all the Ai(s) are obtained immediately and, what is very important, 
always in the most suitable fashion for practical evaluation (as explained before). 



(B.l) 
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Figure 1: The energy as a function of the radius for a specific choice of 



